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ABSTRACT 

It is generally believed that cosmological Gamma Ray Bursts (GRBs) are produced 
by the deceleration of relativistic objects with T > 100. We study the possibility that 
some GRBs are produced along with relativistic matter ejection from supernovae. In 
this model, it is quite likely that the matter has to travel through the progenitor's 
thick envelope before generating GRBs. Under the assumption that the ejected matter 
is described as a single collective matter, we obtain constraints on the matter to have 
r > 100 at the breakout of the progenitor. One advantage of considering this type of 
model is that the expected GRB energy is sufficiently large, in contrast to the GRB 
generation model by the shock breakout in the energetic supernova explosion. We find 
that in general the cross section of the matter has to be very small compared with 
the progenitor's radius and thus the matter has to be bullet (or jet)-like rather than 
shell-like. 



Subject headings: gamma rays: bursts — relativistic shock — stars: supernovae: 
general 



1. Introduction 

It is generally believed that cosmological Gamma Ray Bursts (GRBs) are produced by 
relativistic matter with Lorentz factor F > 100 when the kinetic energy of the matter is converted 
into radiation by inelastic collision with other relativistic matter — internal shock model — , or by 
collision with interstellar matter (ISM) — external shock model (e.g. Piran 1999). 

Although it is still controversial, there are indications that some GRBs might be associated 
with type Ib/c supernovae (SNe) (e.g., Galama et al. 1998, Bloom et al. 1999, Castro-Tirado 
& Gorosabel 1999, Reichart 1999). Especially a type Ic SN, SN1998bw, which was discovered 
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in the error box of GRB980425 has been considered to be the most promising candidate for the 
GRB-SN association (Galama et al. 1998, Iwamoto et al. 1998, Nakamura 1999). In Iwamoto 
et al. (1998) we considered the following scenario for the GRB formation: at the shock breakout 
in the energetic explosion, the surface layer becomes highly relativistic, and the interaction with 
ISM may generate a GRB. One difficulty with this kind of model was that the total energy stored 
in the relativistic region was very small, and hence inconsistent with the observed GRB energy 
(see also Woosley, Eastman & Schmidt 1999). Therefore, something more than the energetic 
explosion is necessary to produce a GRB from a SN explosion. This led to the investigation of 
jet-like explosion (Khokhlov et al. 1999) or the effects of disk accretion around a black hole (BH) 
(MacFadyen & Woosley 1999). 

In this Letter, we assume that during or after core collapse, a small fraction of the matter 
acquires a large energy, becomes relativistic deep inside the SN progenitor, and is ejected. This 
type of matter ejection may be related to the "mystery spot" of SN1987A (e.g., Rees 1987, Piran 
Sz Nakamura 1987, Cen 1999). We show that such relativistic matter can generate sufficiently 
energetic GRBs if the cross section of the matter is sufficiently small. 

The collision of the relativistic matter with the surrounding matter produces an external 
shock, which decelerates the matter. If the relativistic matter is in an optically thin region, 
prompt radiation from it may be observed as a GRB. This kind of situation may be realized, 
if the collapsed star is rapidly rotating and forms an accreting BH, and if matter is sufficiently 
evacuated along the rotational axis as in the model of Fryer &: Woosley (1998). However, numerical 
simulations of such system (MacFadyen Sz Woosley 1999) suggest that although the density is 
significantly reduced near the central BH along the polar direction, it is still high. Also the central 
BH and accretion disk are still surrounded by the envelope of the progenitor. Furthermore, light 
curve modeling of SN1998bw suggests that the SN explosion was globally well described by the 
usual SN model, i.e., the explosive SN shock should propagate through the normal stellar envelope 
(Iwamoto et al. 1998, Woosley et al. 1999). Therefore, it is likely that the relativistic object 
produced near the center has to travel through the optically thick envelope of the progenitor. The 
key question is, thus, if the object can still have F ~ 100, when it breaks out of the envelope. By 
considering a simple model, we obtain the constraints for a relativistic object produced in the 
deep interior of a SN progenitor to have F ~ 100 when it breaks out of the envelope. 



2. Basic equations 

Here, the basic equations and some solutions are shown to describe the slowing down of a 
relativistic object via collision with the surrounding matter. 

We assume that the slow-down of the relativistic object occurs by a series of infinitesimal 
inelastic collisions between the object and infinitesimal masses. We denote the rest-frame energy 
(rest-mass and thermal energy) of the object located at radius r by M{r), its Lorentz factor by 
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r(r), the mass of the envelope that has aheady coUided with the object by m{r), and the thermal 
energy produced in the collision by E. Then the energy and momentum conservation yield (e.g., 
Piran 1999), 

and dE = {T- l)c^dm. (1) 



r2 - 1 M 

The variation of M is given by dM = (1 — e)dE/c^ + dm, where e is the fraction of the 
shock-generated thermal energy that is radiated. When the object is in the envelope, we assume 
e ~ 0. If Tj ^ 1 and F » 1 these equations can be easily integrated to give: 

M(r)^Mo(^)'"'\ (2) 
where Mq and Tj are initial rest-mass and Lorentz factor, and 

m(r)c.^(r(r)-2+^-rr2+^). (3) 

After breaking out from the envelope, the object emits radiation as it is slowed down. If the 
collisions between relativistic objects are not efficient, then slow-down occurs mostly via collision 
with ISM (external shock model). Eq.(^) shows that, by the time the object lose half of its 
kinetic energy (F ~ Fj/2), the mass of m^ ~ m{R)/T has been plowed. The total thermal energy 
produced is (Meszaros & Rees 1992), 

E^ / Tc^dm ^ m{R)c^, (4) 
Jo 

which is assumed to be roughly equal to the total GRB energy. 



In the external shock model, the characteristic synchrotron energy is given by hv = 160 
keV €]i^el(T/100)^n^/^, where eg is the ratio of the magnetic field energy density to the total 
thermal energy, Cg is the fraction of the total thermal energy which goes into random motions of 
the electrons, and n is the ISM number density per cm^ (e.g. Piran 1999). Therefore, F > 100 is 
required for efficient gamma ray emission. 



3. Conditions for GRB generation 

We assume that the objects become relativistic near the collapsed core, and travel through the 
envelope of the SN progenitor. In the following, we consider two cases for the deceleration of the 
objects. One is when the relativistic objects do not lose their energy significantly in the envelope 
of the progenitor, i.e., Fj ~ Tj, where Fj is the Lorentz factor at breakout. The second case is 
when the objects are initially highly relativistic and slow down to Fj <C Fj through the envelope. 
In this paper, for simplicity, we only consider the case with Fj ~ 100. In both cases we find that 
the cross sections of the relativistic objects must be very small compared with progenitor's radius. 
Hence the object should be bullet (or beam or jet)-like, rather than shell-like. Note that Eq.(l) 
applies only when the object can be treated as collective matter, and hence it does not apply to 
long beams or jets. 
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3.1. Case I: Ti ~ T/ ~ 100 

The energy E required to accelerate N' pieces of bullets with mass Mq up to Tj is 
E = N'Moc^Ti. Therefore, 

N'Mo = 5.5 X lO"^^5irri(M0), (5) 

where E^i = £'/(10^^ erg). As shown in the previous section, when each bullet plows a mass of 
~ Mq/T, r becomes roughly half. Therefore, if the bullet is not to be significantly decelerated 
after traversing a mass Menv and radius R, the diameter of the bullets db (here the diameter is 
assumed to be constant) should satisfy the condition: db < {AirMo/SMenv^iy^'^R- For example, 
for the progenitor model of SN1998bw, M^nv ~ IOM0 (Iwamoto et al. 1998), using Eq.(|), we 
obtain 

db < io-y£;5i/iV'(ioo/r)i?. (6) 

The situation considered here is similar to the model considered in Heinz &; Begelman (1999). 
They assume that the central engine of the burst releases A'^' bullets distributed over an opening 
angle of ~ 10° with Fj ~ 1000. Each bullet is assumed to be expanding sideways with a velocity 
Vs = ac/T <C cT measured in the observer's frame (ac is the expansion velocity in the comoving 
frame). Since the viewing ang le 1/r < e, the observed number of bullets is ~ N'/{e'^T'^). They 
considered this model to explain the short-time variability of the canonical GRBs, and claimed 
that the external shock model can explain that variability if A^ ~ 100, a ~ 0.01, F ~ 10^, and 
the ambient gas density n ~ lO^cm"^. In this paper we do not attempt to produce short time 
variability, since it is not certain that a GRB associated with a SN should show it. 

Each bullet may expand if it has sufficiently large internal energy, though the expansion 
rate depends on several assumptions. If we assume that the bullets expand sideways with 
comoving velocity ac, as in Heinz & Begelman (1999), then condition (^ is somewhat relaxed as 
shown in the following. The cross section of the bullet in the observer's frame S{r) is given by 
S{r) ~ (2fstF^)^ = (ar/F)^, where Vg = ac/T is the sideways velocity in the observed frame, and 
t ~ r/(2cF^) is the observer's time when the bullet is located at radius r. The total mass plowed 
by the bullet nis is given by, 

rus = I p{r)S{r)dr, (7) 

where R is the radius of the progenitor, ro(^ R) is the radius where the bullet is accelerated to 
F = Fj, and p{r) is the density of the progenitor. 

If p oc r'^{h < 3), and tq < R, Eq.(0) yields tUs ~ {^71^^ Menva'^T''^ , where Menv 
= J^4:7rr^p{r)dr is the progenitor mass above radius vq. Note that rUg is independent of b. A 
constraint on a is given by nis < Mq/T: 
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The diameter of the bullet at breakout {r = R), d{R) is 

d{R) = aR/r < 10-\ -—^^i——R. 

After the bullet breaks out of the envelope, a GRB is produced. If it is produced by 
the external shock, total energy of the order of Egrb ~ M(? is released as the bullet plows 
through a mass of ~ M/Tf. Here M is the rest-energy of the bullet at r = i2, which is 
estimated by Eq.(|), M ~ MoTi/Tf ~ Mq. Using Eq.(|), we obtain the total GRB energy 
Egrb ~ IQ^'^ E^i{im/T f){N /N') erg. Note that if many bullets {N' » 1) are ejected within 
the opening angle 9, the total GRB energy seen by an observer is reduced by the factor of 
N/N' ~ l/{eTf. 



3.2. Case II: » T/ ~ 100 

Suppose again that the bullets expand sideways with constant velocity a measured in the 
comoving frame. Then the diameter of the bullet at radius r in the observed frame is given by, 

d{r) = r a/T{r')dr'. (10) 



In general T{r) is a function of a and p{r). However, here we simply assume that T{r) = rj(r/ro)~" 
and p{r) = po{r/r(j)~'' in order to get a rough constraint without obtaining an exact solution for 
r(r, a, p{r)). Here the subscript represents the quantity at r = rg in the observer's frame. Strictly 
speaking, assuming that both a and h are constants is not self-consistent, however the result given 
in Eq.(|T^) is not so sensitive to the value of a and h. Therefore, this assumption is sufficient for 
a qualitative analysis. Now Eq. ( |lO|) gives d oc r^"*"". Hence we denote d{r) = dQ{r /rof''^"' , where 
do = otrQ/[{a + l)rj]. If constants a and h satisfy < & < 3, and 2a + 3 — 6 > 0, then Eqs.@ and 
(0) give 

rUs ~ r — — T- (11) 

We define the envelope mass by M^nv = ^j^^TTr"^ p{r')dr' ~ {AtttqPq/^ — b){R/ro)^~^. Now, 
the diameter of the bullet d in the observed frame at the progenitor surface (r = R) is given by 



. iRV^" l 2a + 3-b27rMoTi OMR 2a + 3-b E51 

d[R)-do{^^J -^^y Me„. ^ V 3-6 N'M^UMeY ^ ' 



or 



" = + 'V 3-6 iV'M_(Me)- ^^^^ 

Note that if a is significantly smaller than the right hand side of this equation. Case I applies. This 
result is not sensitive to the values of a and 6 as long as 2a -1-3 — 6/3 — 6~ 0(1)- This justifies our 
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approximate analysis assuming that a and b are constants. The total GRB energy (in the external 
shock model) can be estimated as Egrb ~ W^^ E^iiWO/T f){N/N') erg, which is same as Case I. 

4. Discussion 

We have investigated relativistic matter ejection from SN progenitors in order to study 
the possibility that SNe generate GRBs. If we consider this type of model, it is likely that the 
relativistic matter has to travel through the optically thick envelope of the SN progenitor before 
generating the GRB. We have obtained conditions for such relativistic matter to have T > 100 
after passing through the thick envelope of the progenitor. In general, in order to satisfy these 
conditions, the cross section of the objects must be small. Typically the diameter of the objects 
must be at least factor of W^Ec^^^'^ smaller than the radius of the envelope when they break out 
of the envelope. Therefore, the objects should be bullet (or beam or jet)-like rather than shell-like. 
Here we should note that the opening angle of GRB can be significantly larger than the size of 
each bullet at breakout, because the size of the bullet may increase significantly by the time of 
GRB emission owing to its expansion or to merging with other bullets. 

We make no attempt to explain how such bullets may be produced during core collapse 
since the explosion mechanism of massive stars are not well understood. Indeed, currently no 
mechanism is known to make such bullets. It is only speculation but such bullets may be produced 
if the exploding core is inhomogeneous (e.g. Burrows & Hayes 1996) and the explosion energy is 
concentrated in some small regions, or if the collapsing core forms an accreting black hole and 
relativistic bullets or jets are produced around the black hole (a jet of short duration may look 
like a bullet). 

Next let us estimate the initial rest-frame density of the bullets. Suppose the bullet 
has a volume of roughly in the observed frame, and has T = Fj at mass coordinate 
Mr = 2Mq in a SN progenitor. If we take a progenitor model of SN1998bw, at 
Mr = 2Mq the radial coordinate is rg = 2 x 10® cm. From Eqs. (P) and (|l^), we find 
that do ~ W-^T-^{E5i/N'Menv{MQ)y/'^ro if we require F/ > 100. Here F^ > 100 for Case I 
and Fq = Fj > 100 for Case II. The bullet mass is Mq ~ Fjdgpj. Thus, using Eq.(|5|), one obtains 
a constraint pi > W^Ta{N' /E^iy^^Melv{MQ){gcm^'^). This density is much larger than the 
progenitor's density at this radius /o(ro) ~ lO^g cm~'^, which means that the bullets must be 
produced deeper inside. 

The GRB energy produced by each bullet is Egrb ~ W^^E5iN'~'^{100/T j) erg (in the 
external shock model). If E^i < 1, this is sufficiently large and not too large to explain the energy 
of GRB980425 (~ 10^® /& erg, if its distance is ~ 40Mpc, where < 1 is the beaming factor) 
which might be associated with SN1998bw. 

The GRB980425 lasted about 30 seconds and this must be explained in our model, if it is 
produced by bullets ejected from the progenitor. The time scale of the GRB in this model is 
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determined by the time scale of the deceleration of the bullets, or by the life time of the central 
engine. If it is determined by the latter, we can assume anything because almost nothing is known 
about the central engine. In the following, assuming the external shock model, we discuss how the 
time scale is determined in the former case. 

In our model, the bullet mass at breakout is M ~ MoTi/Tf. Each bullet emits a GRB as it 
plows ISM mass of nig ~ M/Tf. The mass plowed depends on ISM density and on the bullet's 
cross section. If the ISM is produced by wind mass loss from the progenitor, the number density 
can be estimated by nisMii") = 1-5 x 10^^M_6^^20^'"^^ cm~^, where M_6 is the mass loss rate in 
units of lO^^MQt/r^^ and V20 is the wind velocity in units of 20km s^"*^. M-q/v2o is roughly ~ 0.1 
for O stars and ~ 1 — 100 for red giant stars. The time evolution of the bullet size is uncertain. 
The radius may increase due to internal expansion or by merging with other bullets. Here we 
assume that the diameter of a bullet increases as oc r^^^ with k ~ 0(1). Suppose the diameter of 
the bullet is ~ 10~^(-E'5i/A^')^/^rj^i? at r = i?, where R is the radius of the order of the progenitor 
radius. For the model of SN1998bw progenitor, R ^ lO^^cm. Then the time scale of the GRB 
duration, r, is estimated by 

AR Rn ( 10^1 \^ 

sec, (14) 



2T)c T) V^-6^2"o'^ii/ 

where Eqs. ^ and ^ are used, Ai? is the distance the bullet travels during GRB emission, 
and = i2/10^^cm. This result shows that if A; = 0, r is independent of R and typically too 
long: r > 10^(100/rj)^ sec. On the other hand, if A; ~ 0.5, r ~ 30 sec is possible. If k is larger, 
even a milli-second time scale may be possible. Therefore, we may conclude that GRB durations 
of milli-seconds to months may be realized with this model, depending on the parameters and 
situations. Especially, the expansion rate of the bullet radius is very important to determine the 
time scale. In this Letter the expansion rate was treated as a free parameter, and constrained 
it for successful GRB generation. However, the rate may be estimated if the pressure inside and 
outside the bullets is given. The evolution of bullets including their expansion rate will be further 
explored by performing hydrodynamical simulations in future work. 

Finally we comment on the internal shock model for the GRB generation. In the last 
paragraph, we mentioned that the cross section of the bullets may be increased via merging with 
other bullets. If this is the case, the GRB may be produced by the internal shocks produced 
during merging. If the merging bullets grow larger than the angular diameter of F"^ before losing 
their kinetic energy, their subsequent evolution will be similar to the model by Kumar &: Piran 
(1999). They considered collision of relativistic blobs with total kinetic energy 10^^ erg as the 
source of GRBs, and showed that the observed diversity of the GRB energy can be explained by 
the inhomogeneity of the blobs and the differences of lines of sight. In general, the efficiency of 
the kinetic to thermal energy conversion is much worse in the internal shock model than in the 
external shock model. However, the time scale of energy loss can be several orders of magnitudes 
smaller than the external shock model and this may be an advantage (e.g. Piran 1999). 
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